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Abstract It is shown that time reversibility of Hamiltonian microscopic dy- 
namics and Gibbs canonical statistical ensemble of initial conditions for it 
together produce an exact virial expansion for probability distribution of 
path of molecular Brownian particle in a fluid. This expansion leads to in- 
equality connecting logarithmic derivative of the distribution with respect to 
fluid density and characteristic volume occupied by pair statistical correla- 
tion between the path and fluid molecules. Due to the inequality, finiteness 
of this volume means that asymptotic of the distribution is essentially non- 
Gaussian. For principal case when fluid is dilute gas it is demonstrated that 
the pair correlation volume is actually finite and bounded above. Therefore 
even under the Boltzmann-Grad limit the path distribution possess power- 
law long tails (cut off at distances of ballistic flight). 
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1 Introduction 

Thermal motion of molecules in fluids or quasi-particles in solids gives us ex- 
amples of chaotic walk of small molecular-size "Brownian particle" . As a rule, 
considerations of such walk, or "molecular Brownian motion" , involve more 
or less amount of stochasticity constantly brought in from random "thermal 
bath" or random choices like coin tossing, etc. Many understand that this 
way of thinking conflicts with principles of statistical mechanics where the 
only random choice is that of initial conditions for Hamilton equations. Some 
understand that even the Boltzmann's "Stosszahlansatz" and thus validity of 
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the Boltzmann equation, even for dilute gas under the Boltzmann-Grad limit, 
remain unproved [HE]- All the more stochastic probability-theoretical mod- 
els of molecular motion in fact stay based on faith rather than Hamiltonian 
dynamics. 

Such state of affairs was predicted by Krylov in [3] (firstly this book was 
published in 1950 in Russian). According to him, actual (cause-and-effect) 
independence of events in reality, on concrete phase trajectory of a system un- 
der interest, generally does not imply statistical independence of same events 
in theory, in statistical ensemble of trajectories induced by a given ensemble 
of their initial conditions. For instance, ensemble-average density of pair colli- 
sions of gas atoms, along with pair distribution function of colliding atoms, in 
spatially non-uniform gas does not reduce, even under the Boltzmann-Grad 
limit, to (quadratic form of) one-particle distribution function [4]. 

In such case kinetic theory loses its traditionally most powerful tool. To 
compensate this loss, one can examine other tools, e.g. exact "generalized 
fluctuation-dissipation relations" (FDR) [6ll7l[8] implied by determinism and 
time reversibility of microscopic dynamics. In the present paper we exploit 
one of simplest FDR what correspond to the classical Gibbs canonical ther- 
modynamically equilibrium ensemble. 

Main object of our interest is probability distribution, Vo{t, AR) , of dis- 
placement, Z\R , of molecular Brownian particle (BP) during time interval 
(0, t) much longer than its mean free-flight time. We transform the FDR into 
an original virial expansion of Vo{t, AR) connecting its n-order derivatives 
with respect to density of gas (or, generally, fluid) and {n + l)-particle prob- 
ability distribution functions and correlation functions describing statistical 
correlations between the BP's displacement and gas particles. As combined 
with positiveness of distribution functions, the virial expansion produces dif- 
ferential inequalities to be satisfied by VQ{t, AR) . We show that already 
first of them imposes a limitation on possible steepness of Vo{t, ARYs tails 
at |R| > V 6Dt , where D is BP's diffusivity, certainly forbidding exponen- 
tially vanishing tails. In particular, the Gaussian asymptotic is forbidden, 
first of all, for the Boltzmann-Grad gas. 



2 Fluctuation-dissipation relations 

Let us consider a molecular-size "Brownian particle" (BP) among iV ^ 1 
atoms in a box with volume f2 under the thermodynamical limit: N oo , 
n ^ oo , N/fl = Vf) = const. Introduce designations q = {R, ri, ... , tat} 
and p — {P, pi, ... , Pat} for canonical coordinates and momenta of our 
system, respectively, and 
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for its Hamiltonian in absence of external fields, so that H{q, — p) — H{q, p) . 
Besides, let q(t) = q(t, q, p) and p{t) = p{t, q, p) be solution to correspond- 
ing Hamilton equations under initial conditions q(0) = q and p(0) = p , and 
/9eq(q, p) cx cxp [— iJ(q, p)/T] equilibrium probability distribution for these 
conditions (normalizing factor is omitted). Then for any suitable functions 
A(q, p) and B{q,p) we can write 

A(q(t),p(t))B(q,p) peq{q,p)dqdp = (1) 




B{ci{t),~p{t)) A(q, -p) pe?(q, P) dqdp 



This quite obvious identity is very particular form of the mentioned FDR. 
But we want to concretize it even greater by choosing 

A(q,p) = <5(R-R') , 

B(q,p) = /2<5(R-Ro)<5(P-Po)exp[-E^.C/(rj,p,)/T] (2) 
- 12 <5(R - Ro) <5(P - Po) n, [ 1 + P,) ] , 

where ip{r,p) = exp [— L/(r,p)/T] — 1 . 

We assume that the potentials of atom-atom interaction, Uaa{f) , and 
BP-atom interaction, Uab{f) , satisfy conventional standards of kinetic the- 
ory of fluids [HlQllTO] which allow to speak about solutions of the Hamilton 
equations at arbitrary large N at least in statistical language, in terms of 
the distribution functions (DF) as introduced by Bogolyubov [9| . What is for 
the function U{r, p) , or ip{r, p) , it is good enough if 

J\4>{v)\dr < oo , 0(r) = /^(r,p)G„(p)dp , (3) 

with Gm{p) — (27rTTO)^^/^ exp (— p^/2Tto) being equilibrium Maxwell 
momentum distribution. At these assumptions the thermodynamical limit 
of H]), when the box walls move away to infinity, can be treated in full anal- 
ogy with J9]. Firstly, consider right-hand side of ([T]). 




3 Many-particle distribution functions 
and correlation dressing of BP 

The right-hand side of ([l]), after substitution of tends to 

J J B{q{t),-p{t))A{q,-p)peq{q,p)dqdp - Rq, -Pq, |R'} , 

.F{i,Ro,Po, V^IR'} = Fo(t,Ro,Po|R') + (4) 
+ / ^n(^,Ro,ri...r„,Po,Pi...p„|R')n^(^J'-P:') ' 

n=l "''■XP j=l 

where J"^ means integration over n coordinates ri...r„ and momenta 

pi...p„ ; function Vb(<, Ro, Po|R') is conditional probability density of 
finding BP at t > at point Rq with momentum Pq under condition 
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that BP had started at < = from point R' ; function F„ is joint condi- 
tional probabihty density of this event and simultaneously finding some n 
atoms at points Tj with momentums Pj under the same condition. 

In respect to atoms, all these are completely analogous to usual non- 
normalized DF of infinite gas [S]- The role of normalization is played by 
uncoupling of inter-particle correlations at infinity; 

F„(... Tfc ... pfc ... ) ^ F„_i(... rfc_i, Tfc+i... pfe-i, Pfe+i... ) Gm(pfe) , 
Fi(t,Ro,Po,ri,pi|R') -> yo(i,Ro,Po|R')G,n(pi) , ^ ' 

when Tfc — > oo and ri — s- cx) , respectively. But in respect to BP all DF from 
are normalized in literal sense. In particular, 

/ l/o(i,Ro,Po|R')c?Po = Fo(t,Ro-R') , / Fo(i,Ro-R')'^Ro = 1 

In respect to ip ■, expression J-{t, Rq, Pq, V' |R'} represents generating func- 
tional of the DF quite similar to the Bogolyubov's one [9]. Under condition 
(121) the series in ^ converges, hence, this functional is well defined. 
From ^ it is obvious that initial conditions to these DF appear as 

Fo(0,Ro,Po|R') = <5(Ro-R')Gm(Po) , 

K(0,Ro,ri ...r„,Po,pi ...p„|R') = 

- F,l'^'?^(ri...r„|Ro)5(Ro-R')GAf(Po)n;=iG™(p,) , 

.F{0, Ro,Po, ^|R'} = (5(Ro-R')Gm(Po)^^^''H'/'|Ro} , 

n=l 3=1 

where functions Fi'^''^(ri... r„|Ro) are conditional thermodynamically equi- 
librium DF of the fluid, corresponding to condition that BP is placed at point 
Ro , and !F'^^'^\4>\R.q\ is generating functional of these DF. 

Notice, first, that Fi'^'^(ri... r„|Ro) and thus Fn{t — 0, ...) take account 
of all correlations between BP and atoms, as well as atoms themselves, what 
are caused by their interactions. Second, Fn{t > 0, ...) treat essentially 
the same equilibrium state as at t = Q . Nevertheless time dependence of 
Fn{t > 0, ...) does not copy that of Vo(i, Ro, Po|R') , i-e. 

Fn{t, Ro, l"l ■■■ fn, Poi Pi ■■■ Pn |R') 7^ 

^ (i , Ro , Po I R') -F^n^'^i-i . . . r„ I Ro ) n, G„^ (p, ) 

Hence, at t > some extra correlations do arise between BP and atoms. 
Statistical meaning of these correlations clears up if notice that in absence of 
information on start position of BP, R' , all joint DF of BP and atoms would 
be invariable. This shows that in fact extra correlations connect total previous 
bp's displacement, or path, Rq — R' , accumulated during time interval 
(0, i) , and current state of the fluid at BP's surroundings. Therefore one 
may name them "historical correlations" . Anyway, regardless of naming, they 
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should be attributed not to separate phase trajectories but to the statistical 
ensemble of trajectories as the whole. 

Let us introduce functions Rq, ri ... r„, Po, Pi ... Pn|R') which de- 

scribe "historical correlations" in a pure form. By tradition, differences of 
time- varying DF from equilibrium (or quasi-equilibrium) DF are termed "cor- 
relation functions" (CF) ^lll9[|TT]. We will apply this term also to the "his- 
torical correlations" and define corresponding CF through their generating 
functional as follows: 

.F{i,Ro,Po, ^|R'} = T^^^Hmo} V{t,-Ro,PQ, ^P\-R'} , 

V{i,Ro,Po, V'lR'} - Vo(i,Ro,Po|R') + (6) 

n=l "''■XP 3=1 

with (j) = (j)(r) introduced in ([3]). At n = 1 , in particular, 

Fi(i,Ro,ri,Po,Pi|R') = (7) 
= K)(i,Ro,Po|R')^i''^(i-i|Ro)G™(pi) + Fi(i,Ro,ri,Po,Pi|R') , 

where function Vi describes pair historical correlation of BP with atoms. 

According to definition ^ and above written initial conditions for DF, 
all the CF vanish at start of BP's observation: Vn{t = ,...) = {n > 0). 
Next BP gradually gets dressed in correlations with the fiuid. In view of ([5]) 
and (O, Vn{t , ...Vk...) ^ at Vk oo , that is the correlations always 
vanish at infinity. Below we will see that their actual spatial spread is very 
important characteristics of the Brownian motion under consideration. 



4 Response of BP's path to fluid 
perturbations and its relations with CF 

Further, consider left side of identity ^ under choice ^ rewriting it as 

// A(q(t),p(t))B(q,p)/9eg(q,p)dqdp = / / !f (q, p) Pe^Cq, p) dqdp x 
X / / 5{R{t) - R') <5(P(0) - Po) p»„(q, p) dqdp , 

(8) 

where !Z'(q, p) = S(R - Rq) [ 1 + V'(ri, Pj ) ] and 

Pzn(q,p) = >Z'(q,p)pe9(q,p) [//f(q',p')Peg(q',p')f^q'c^PT^ 

Evidently, last integral in ([S]) represents averaging over new statistical ensem- 
ble of initial conditions, with new probabilistic measure Pm(q, p) . The lat- 
ter is thermodynamically non-equilibrium measure which however would be 
a canonical equilibrium one in presence of imaginary (momenta-dependent) 
external potential U (q, p) = — T In [1 + p) ] . Hence, 

/ / 5{R{t) - R') 6{P{0) - Po) p.„(q, p) dqdp = , . 

= T/{t,R'|V,Ro,Po}GM(Po) , ^ ' 
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where V{t, R'lV', R-Oi Po} is conditional probability density of finding BP at 
t > at point R' under conditions that initially, at t = , it was located at 
point Rq with momentum Pq while the fluid was in such non-equilibrium 
spatially non-uniform state what turns into equilibrium under the external 
potential U{q, p) . Noticing, besides, that 

//if(q,p)pe,(q,p)rfqrfp = ^('"^^0} , 
we can transform equality ([8]) to 

J J A{q{t),p{t)) B{q,p) peqiq,p) dqdp = , , 

= y{t,R'|V,Ro,Po} Gm(Po) ^(^«H0|Ro} ^ ' 

Instead of i^iq, p) or U (q, p) , the non-equilibrium ensemble determined 
by Pm(q, p) can be characterized in terms of corresponding initial mean 
densities of atoms in the /i-space, /i{r, p| ?/;, Rq} , and coordinate space, 
z^{r| 0, Rq} — J /i{r, p| ip, Ro} dp . More precisely, that are conditional den- 
sities at given BP's position. In other words, /J,{r, p| -0, Ro} is initial one- 
particle DP of the fluid. It is easy to verify that 

At{r,p|0,Ro} = [l + VXr,p)] G„,(p) ^ ^'"'^^'^1^^'^°^ (H) 

After all, combining formulas fl]), ([1]), (l6|) and (fTO|) . we come to our main 
formally exact relation: 

y{i,R'|V',Ro,Po} Ga/(Po) = T4,(i,Ro,-Po|R') + (12) 
+ K(i,Ro,ri...r„,-Po,Pi...p„|R') n^^^'j'^Pj) 

n=l "''■XP j=l 

It connects, from left hand, probability distribution of BP's path in initially 
non-equilibrium nonuniform fluid and, from right hand, similar distribution, 
along with generating functional of correlations between the path and current 
bp's environment, for equilibrium uniform fluid. At ip{r,p) = this is 
merely time symmetry relation for equilibrium random walk: 

yo(i,R'|Ro,Po) Gm(Po) = yo(i,Ro,-Po|R') (13) 

Of course, all the DF and CP are dependent on the mean gas density 1^0 . 
But for brevity in relation p2)) . as well as before it and almost everywhere 
below, corresponding argument is not written out. 



5 Virial expansion 

Let us choose ijj{r,p) independent on momentum, tjj{r,p) — (j){r) , with 
0(r) being constant, 0(r) = = const, inside some sphere |r — Ro| < ^ and 
vanishing outside it in a suitable way. Then factual perturbation of spatial 
uniformity of the fluid and thus of its equilibrium initially takes place at 
|r-Ro| > C only. 
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Condition ([3]) requires from the radius ^ to be finite. Nevertheless we can 
make it as large as we want. Therefore, at any given finite time to we can 
find such finite f that at t < to the whole BP's path on left side of for 
sure realizes in uniform equilibrium fluid, although with a density i' which 
differs from the density vq on the right-hand side. For instance, already such 
value as ^ = kvsto, where Vg is speed of sound in our fluid and fc > 2, 
seems to be sufficient. 

Under these assumptions, the functional R'|(/i, Rq, Po} simplifies to 
mere function, and for t < to we can write 

/F{t,R'|0,Ro,Po}GM(Po)dPo = K,(i,R'-Ro; t') , (14) 

where third density argument, i> , of the BP's path distribution is introduced, 
so that Vo{t,Ro - R' ; i^o) = K)(i,Ro - R') • On right-hand side of (dH), 
under same conditions, none correlation between BP and fluid may propagate 
up to distances |r — Ro| > ^ . Therefore in all integrals 0(r) can be replaced 
by the constant. 

These reasonings show that relation tolerates extension to ip — (p = 
const . At that after use of and (fH)) it transforms to 



(h 

Vo(t,Z\R; i^{i^o,<l>)) = Vo{t,AR; i^o) + V ^Vn{t,AR; vo) (15) 

n=l 

Here AR = Rq — R' , functions Vn with n > are defined by 

Vn{t, AR-jVo) — 1^0 / Kt(i, Ro, Ti... r„, Pq, Pi... Pn I R') dPo (16) 

J rxp J 

and v^uo^ip) is the left fluid density v as a, function vo and 0. Evidently, 
that is z^{r|(/),Ro} taken far from BP, at |r — Ro| 3> Vb , where rf, is char- 
acteristic radius of BP-atom interaction potential Uab{r) ■ Due to (fTTjl . 

i^ = Ki^o,0) = i^o{l + cl>){l + iyocf>J[F^"'\r)-l]dr + (17) 
+ '^j\Ft'\r,,r,,0)~Ft''\r,)-F^^^\r,)~Ft'\r,-r,) + 2]+...} 

with F^'"'\r) and F!^'"^\ri , r2,r3) being standard equilibrium pair and 

triple DP at density vo j and symbolizes integration over ri and r2 . 

Formula (fT5|) may be rated as virial expansion of the BP's path prob- 
ability distribution. But it signiflcantly differs from the well-known virial 
expansions of thermodynamical quantities |12| or kinetic coefficients jl3j : 
the latter expand over absolute value of density i/ while (fT5|) in fact over 
difference v/vo — 1 or logarithm ln(z^/i^o) • 

For further, the first-order term of the virial expansion is most important. 
Differentiation of (fTSl) with respect to , at 4> — , together with (fT6| and 
p7)) yields 

^ dVoit,AR; vo) /,R ^ n«^ 

vo ^ ^ Vi{t,AR; vo) = (18) 

ovo 

^ vo / / Fi(t,Ro,Po,r,p|R') dr dPodp , 
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= -0 + I (r) - 1 ] dr = .0 T (^1^) ^ (19) 

The last equality, where P denotes pressure, is known from the classical 
statistical thermodynamics [12] . 



1^0 



6 Restrictions on BP's path distribution 

Relations p5|) with (fTB|) and (fT5|) with (fT5|) prove that historical correlations 
between BP and fluid really exist. Let us show that they imply essential 
restrictions on possible shape of the distribution Vb(t, Z\R \vo) ■ 

Consider the pair correlation Rq, Pq, r, p |R') . Let from now fi 

denotes a finite region in the (r — Ro)-space and at once volume of this 
region. Introduce Q{5) = fi{t, Z\R, Pq, p, 5) as minimum of all such regions 
Q what satisfy 



/,^V^i(t,Ro,r,Po,p|R')dr - / V^i(i, Rq, r, Pq, p|R' ) dr | < 
< 5 |/Fi(t,Ro,r,Po,p|R')cir| 



(20) 



with some fixed < (5 < 1 . Thus Q{5) = ]7(t, ziR, Pq, p, (5) is minimal 
region containing at least 100 (1 — 6) percents of the total pair correlation. 

From the other hand, integrate identity ([7]) over ri e f2{6) . Since (O 
represents a probability density, result must be non-negative: 

Vo{t, Ro, Po|R') G„,(p) F('=')(r|Ro) dr + 

+ Fi(t,Ro,r,Po,p|R')dr > 

This inequality together with (|20p yields more interesting one: 

K)(i,Ro,Po|R')G„(p) /,,(5)F('=«^(r|Ro)dr+ .^l) 
+ (l-<5)/Fi(t,Ro,r,Po,p|R')rfr > 

Next perform here integration over momenta, apply virial relation and 
then divide all by (1 — 6) and multiply by vq . Thus we obtain 

n{t, AR, 5) Vo{t, AR ; uo) + ?o ""^ > (22) 

oiyQ 

Here J7(t, Z\R, S) is characteristic pair correlation volume defined by 

n{t,AIl, S) = (l-(5)-i //GM(i,Po|^R)G„,(p)x 

X Init,ATL.Po.P,5) ^t''\^\^o)dr^ dpdPo 

(23) 

with GA/(t,Po|Z\R) = Vo(i,Ro,Po|R')/^o(i, ^R) being conditional dis- 
tribution of momentum of BP under given value of its path. 

Inequality (|22|) presents us a restriction on possible negative values of 
derivative dlnVo{t, AH; i/o)/9 In j/q . The smaller is f2{t,A'R, S) the stronger 
is this restriction. Therefore n{t, AH, 5) can be replaced by 

Ha, AR) = min TlU, AR, S) 

^ ' 0<<5< 1 ' ' 
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Now let us suppose that at sufficiently large temporal and spatial scales, 
at i ^ r and |Z\R| 3> A, with r and A — vtt being BP's mean free- 
flight time and mean free path, respectively (and vt ^ \/T/M its thermal 
velocity), the BP's path distribution tends to the Gaussian distribution: 

Vo{t,A'R) VG{t,AR) = (ATTDt)-^/^ cxp {- AR^ / 4Dt) (24) 

Then inequality (P^ . with use of (|19p. takes the form 



1^0 nit, AR) + 



dD /Z\R2 3 



VGit,AR) > (25) 



In natural situation when BP's diffusivity is a decreasing function of the den- 
sity, dD/duQ < (e.g. when our fluid is a gas), the second addend in square 
bracket becomes negative at z = AR^/ADt > 3/2. As the consequence, 
inequality (|25|) can be satisfied if and only if at large values of z the volume 
f2(t, AR) infinitely grows, at least proportionally to z . 

Thus, we came to intriguing paradox or theorem. If asymptotical statis- 
tics of the BP's path is Gaussian then characteristic volume occupied by 
pair correlation of BP with atoms is unbounded above. In opposite, if the 
correlation volume is bounded above, i7(<, Z\R) < fimax with f^max being 
independent on t and AR, then asymptotic of Vb(i, Z\R) must be essen- 
tially non-Gaussian (see below). 

What is closer to reality? Further, we want to show that the second 
statement conforms to truth at least if our fluid is a dilute gas. 



7 Pair correlation volume in dilute gas 
and the Boltzmann-Grad limit 

For our purpose, to exclude non-principal complications, the Boltzmann- 
Grad limit (BGL) is very useful, when ^ oo while r;, — > [va 
is radius of short-range repulsive atom-atom interaction) in such way that 
gas non-ideality parameters 47rrj^i'o/3 and 47rr^i/o/3 vanish but mean free 
paths of BP, A = (7rr^^'o)~"'^ , and atoms, A — (Trr^j/o)^^ , stay fixed. 

It is easy to justify that under BGL expression (fTT|) effectively simplifies: 

/i{r,p|i/',Ro}/i/o ^ [l-hV(r,p)]Gm(p) , function F^''''^ (r|Ro) in ^ can 
be replaced by unit, and inequality ((22l) yields 

Q{t,AR) Vn{t,AR; v^) + ^ /^"(^'^^i > q , (26) 

-n{t,AR)= I I GM{t,Vo\AR)GM min ^^i^l^^i^^lM dprfp^ 

J J 0< i5< 1 1 — d 

Let us discuss the pair correlation volume, basing on the first-order term 
of the expansion (fT2|) : 

= i.o"t^i(t,Ro,r,-Po,-p|R') (27) 

■0 = 



Gm(Po) 



(5V-{t,R'|^,Ro,Po} 
(5iA(r,p) 
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At that, we will take in mind also a few things about pair correlations and 
pair CF known from classical theory of dilute gases [T|[9jfT0lfTT] and, besides, 
those distinctive property of the Boltzmann-Grad gas that any two of its 
particles deal with different non-intercrossing collections of other particles, 
except cases when direct encounter of the two takes place. 

Notice that left side of ([27|) reflects reaction of probability distribution 
of bp's path, R' — Rq , to lodging, at t = , in point r one extra atom 
with momentum p . In view of the stated property, this reaction is nonzero 
at such initial relative disposition of BP and extra atom only which ensures 
their direct collision, either soon forthcoming or just happening one. For 
that, two conditions should be satisfied. Firstly, vector r — Rq belongs to 
so-called "collision cylinder" which has radius ~ rj, and is oriented along 
relative velocity of the particles, v — Vq = p/m — Pq/M . Secondly, |r — 
Ro I is not much greater than yl or A . Otherwise collisions of either BP or 
given extra atom with other atoms (the rest of gas) will prevent the desired 
collision. These remarks, in the framework of (|27p . make it obvious that 
volume occupied by the pair correlation can be estimated as 



(assuming for simplicity that X ^ A), i.e. it is comparable with volume 
displayed per one atom of the gas. 

The first condition agrees with the classical theory: the non-zero reaction 
on the left in (jST]) corresponds to non-zero correlation on the right well- 
known at least for post-collision dispositions. However, the second condition 
does not agree: a spread of this correlation along the collision cylinder in 
classical theory remains indefinite. Formally, infinite. The matter is that this 
theory neglects destruction of pair correlation by collisions of the pair with 
other particles. This leads to infinite volume of pair correlation, which is 
highly doubtful result. 

At the same time, the theory [l][9lll01[TT| says that magnitude of pair CF, 
Vl(^, Ro, r, Pq, p|R') , inside colhsion cylinder generally is comparable with 
product of one-particle DF (first right-hand term in ([7])). This is very impor- 
tant statement, and it looks, at least for |r — Ro| < yl, as quite doubtless 
result. It follows also from the first of the BBGKY equations under remark 
that first term of ([7]) does not contribute to the "collision integral" (see [U 
[51 [H[|15j ). In view of the exact virial relation psp this result clearly cancels 
the previous one. Moreover, confirms that in fact volume of pair correlation 
has an order of ■ 

Importantly, this estimate manifests no visible dependence on either t 
and Z\R or the momenta. Therefore, taking into account that at t 3> r 
certainly Gjif (t, Po|/iR) « G'Af(Po) (see [4 l [5 | ll4 p i5j). we can write 



n{t, Z\R,Po,p,<5) 




1 



■o 




(28) 



with f2{S) representing i7(t, ziR, Pq, p, 5) at typical thermal momenta (see 
its definition in previous section). As the consequence, inequality (|26p can be 
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replaced by 



c,Vo{t,AK; .0) + > , (29) 



where ci — ^max ^ 1 is a constant of order of unit. 



8 Probability distribution of molecular 
random walk in dilute gas has long tails 

The Gaussian asymptotic (IM)) is incompatible with inequaUty To see 

what is allowed instead of ((24|) . it is natural to suppose that at i ^ r 
and |^R| ^ A the BP's path distribution Vb(i, Z\R; J^o) is characterized, 
similarly to Vait, AR) , by a single parameter, that is diffusivity. Then 

Vo(t, AR ; uq) (4L't)-3/2 ^{AB?/4.Dt) (30) 

Conditions / dAR = 1 and / AH? dAR = 6Dt (which in fact 
serves as quantitative definition of the diffusivity) imply that 

/ If (a2) da = 1 , / a^if (a2) da = 3/2 

Since in a dilute (Boltzmann-Grad) gas diffusivity is inversely proportional 
to gas density, D — vtA cx Vq"^ , that is Vq dD/dvo ~ , inequality (^5]) 
yields 

. . d^(z) 3 , , 

a^{z) + z > , a = - + Cl (31) 

Hence, function ^{z) must have power-law long tail: ^(z) oc z~ ^ at 
z — > 00 , where /? < a . At that, the second condition requires a > 5/2 . 

One may rightly object that such asymptotic is nonphysical since has 
infinite higher-order statistical moments of AR and makes probable enor- 
mously fast displacement of BP. But we can easy correct all this if recall that 
in addition to diffusivity the BP's motion has one more significant parame- 
ter, namely, BP's thermal velocity vt , which is independent on t^o and not 
touched upon by ((29|l . Hence, a true asymptotic has the form 



AR^ \AR\ 



ADt vrt 



iziRI 



VTt 

(32) 

Here function 0(y) « 1 at y <C 1 and fast enough tends to zero at y > 1 
thus sharply cutting long tails of Vb(t, AR) when |Z\R| > vxt . 
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9 Comparison with earlier results 

To get more details of Vq {t, ATI) 's behavior one has to supplement virial 
relations by additional tools, in principle, by corresponding Bogolyubov- 
Born-Green-Kirkwood-Yvon (BBGKY) equations. Anyway, remarkably, our 
present results are in full qualitative agreement with what was obtained in 
[4] (or see [5]) and in [M] from approximative analysis of BBGKY equations 
for the Boltzmann-Grad gas (with marked gas atom in the role of BP). An 
expression for the BP's path distribution, Vo{t, AH), found in^ has just the 
form ^ with !f (z) = r(7/2) (1 + z)-'^^'^ , that is a = 7/2 , which corre- 
sponds to ci = 1^0 f2 max = 2 . At that, the cut-off function starts to work 
already from fourth-order statistical moment: 

(Z\R4(t)) = J AR-^VQit,AR)dATl ^ 3(ATl^{t))^ In^ (33) 

with (Z\R^(i)) = 6Dt . This agrees with results of [4j found in different 
approach when only four lowest moments were under consideration. 

Such behavior of the fourth-order statistical moment as presented by ([33]) 
can be phenomenologically interpreted as manifestation of scaleless 1 /f fluc- 
tuations in diffusivity (and mobility) of molecular Brownian particle. Thus, 
the present work sustains the theory of microscopic Brownian motion and 
related 1/f-noise in (infinitely) many-particle systems which was first sug- 
gested in works p^lflT] . Heuristic explanations of statistical nature of these 
fluctuations can be found also in [illS HTilfTSlfTO] . 



10 Conclusion 

We derived exact relations which connect probability distribution of path of 
molecular Brownian particle (BP) in a perturbed fluid and specific "histori- 
cal" statistical correlations between the BP's path and molecules of equilib- 
rium fluid. The simplest of this relations produced well rigorous differential 
inequality which says that Gaussian asymptotic of the path distribution is 
incompatible with finiteness of spatial spread of the pair historical correla- 
tion, even in dilute gas under the Boltzmann-Grad limit. Then we argued 
that indeed the spread is bounded above. Consequently, the path distribu- 
tion has essentially non-Gaussian form, with power-law long tails (lasting up 
to |Z\R|/t or order of BP's thermal velocity). 

These results mean that BP's path, even in the Boltzmann-Grad gas 
(moreover, first of all in it), can not be disjointed into some parts what 
are statistically independent in the sense of the probability theory. In other 
words, real microscopic dynamics of (infinitely) many-particle system pro- 
duces much more rich randomness than poor probability-theoretical "dice 
tossing" can do. That is occasion to recall the A. Einstein's words: "God 
does not play dice" . 

I acknowledge my colleagues Dr. Yu. Medvedev and Dr. I. Krasnyuk for 
interesting and useful discussions. 
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